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Robust  Stability  Under  Mixed  Time  Varying,  Time 
Invariant  and  Parametric  Uncertainty 

Fernando  Paganini  * 


Abstract 

Robustness  analysis  is  considered  for  systems  with  structured  uncertainty  involving  a  com¬ 
bination  of  linear  time-invariant  and  linear  time-varying  perturbations,  and  parametric  uncer¬ 
tainty.  A  necessary  and  sufficient  condition  for  robust  stability  in  terms  of  the  structured  singular 
value  pi  is  obtained,  based  on  a  finite  augmentation  of  the  original  problem.  The  augmentation 
corresponds  to  considering  the  system  at  a  fixed  number  of  frequencies.  Sufficient  conditions 
based  on  scaled  small-gain  are  also  considered  and  characterized. 


1  Introduction 

A  substantial  amount  of  research  in  recent  years  has  been  devoted  to  analysis  and  synthesis  of 
control  systems  to  achieve  robust  stability  and  performance  in  the  presence  of  structured  uncer¬ 
tainty.  This  implies  a  decentralized  nature  of  the  uncertain  perturbation,  which  is  a  reasonable 
modeling  choice  for  complex  systems,  where  uncertainty  may  be  introduced  at  the  subsystem  level 
(see  Safonov  [17]  and  Doyle  [5]  for  early  treatments  of  this). 

In  addition  to  this  “spatial”  structure,  different  assumptions  can  be  made  on  the  dynamic 
properties  of  the  uncertainty:  real  parametric,  linear  time  invariant  (LTI),  linear  time  varying  (LTV) 
or  nonlinear  perturbations.  All  these  uncertainty  classes  arise  naturally  in  modeling.  Parametric 
uncertainty  appears  frequently  in  first  principles  models;  LTI  perturbations  are  well  suited  when 
‘Electrical  Engineering,  M/S  116-81,  California  Institute  of  Technology,  Pasadena,  CA  91125 
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there  is  frequency  domain  information  about  the  system,  to  be  “covered”  by  a  suitable  ball  of  LTI 
uncertainty.  Time-variations  can  be  captured  by  LTV  uncertainty;  more  generally,  an  arbitrary 
LTV  operator  is  equivalent  to  a  norm  constraint  between  signals,  and  therefore  provides  very  crude 
uncertainty  information  which  may  be  used  to  cover  a  contractive  nonlinearity  in  the  system. 

In  [5],  robustness  analysis  for  structured  uncertainty  was  focused  on  constant,  complex  matrices, 
and  led  to  conditions  in  terms  of  the  structured  singular  value  fi\  an  evaluation  of  (complex)  fx  across 
frequency  captures  the  case  of  LTI  perturbations  [11].  Subsequent  work  [7,  19]  has  considered  a 
combination  of  uncertain  parameters  and  LTI  perturbations,  and  leads  to  mixed  (real/complex) 
^-analysis,  These  conditions  are  exact,  although  their  evaluation  is  computationally  hard  (see 
[4])  and  is  therefore  usually  approached  by  means  of  bounds.  In  particular,  a  class  of  sufficient 
conditions  for  robust  stability  can  be  stated  in  terms  of  a  scaled  small-gain  theorem,  which  can  be 
evaluated  via  convex  optimization.  The  scales  are  chosen  to  commute  with  the  spatial  structure  of 
the  uncertainty,  and  can  either  be  constant  or  varying  in  frequency. 

Robustness  analysis  under  LTV  uncertainty  has  seen  major  progress  in  recent  years,  since  it 
was  discovered  that  the  above  mentioned  convex  tests  apply  exactly  to  this  case.  It  was  found  in 
[8]  for  the  foo  setting,  and  later  in  [18],  [9]  for  the  l^  case,  that  the  constant  scales  condition  is 
necessary  and  sufficient  for  robust  stability  under  structured  LTV  perturbations.  In  reference  to 
the  frequency-varying  scales  test,  it  was  shown  in  [14]  (see  also  [15])  that  this  condition  is  necessary 
and  sufficient  for  robust  stability  against  arbitrarily  slowly  varying  (a  slightly  larger  class  than  LTI) 
uncertainty. 

This  paper  considers  the  situation  where  a  mixed  structure  of  LTV,  LTI  and  parametric  per¬ 
turbations  affects  the  system.  This  is  a  very  natural  problem  since  the  same  reasons  which  lead 
to  “decentralized”  uncertainty  will  often  produce  uncertain  models  involving  a  combination  of  the 
above  mentioned  classes. 

The  inclusion  of  time-varying  uncertainty  precludes  a  straightforward  frequency  domain  fi- 
analysis  as  in  mixed  LTI/parametric  problems.  For  this  reason,  we  propose  an  augmentation 
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procedure  across  a  number  of  frequencies,  which  is  related  to  a  lifting  technique  for  fi  analysis  by 
Bercovici  et.  al.  [2],  and  to  a  power  distribution  Lemma  from  Poolla  and  Tikku  [14].  The  main 
result  in  this  paper  is  to  show  that  a  fi  test  in  the  augmented  structure  is  necessary  and  sufficient  for 
robust  stability  under  mixed  LTV/LTI/parametric  uncertainties.  The  augmentation  also  provides 
an  alternative  for  the  formulation  of  convex  upper  bounds,  and  leads  to  an  exact  characterization 
of  these  tests. 

The  paper  is  organized  as  follows.  In  Section  2  we  review  in  more  detail  the  background 
material  described  above.  In  Sections  3  we  focus  on  the  combined  LTV/LTI  problem,  and  obtain 
the  complex  /i  test  for  analysis.  In  Section  4  we  analyze  the  properties  of  the  corresponding  convex 
upper  bounds.  In  Section  5  parametric  uncertainty  is  also  introduced  in  the  problem,  and  an 
augmented  mixed-^  test  is  derived.  Section  6  contains  examples  which  demonstrate  the  results, 
and  conclusions  are  given  in  Section  7.  Preliminary  versions  of  these  results  were  presented  in  [13]. 

2  Preliminaries 

Robust  stability  and  performance  analysis  under  structured  uncertainty  has  been  the  focus  of 
a  substantial  research  effort  in  recent  years.  This  section  contains  a  (by  no  means  exhaustive) 
summary  of  previous  work  related  to  this  paper. 

A  standard  setup  for  robustness  analysis  is  depicted  in  Figure  1.  This  picture  represents  a 
robust  stability  problem,  where  M  is  the  nominal  system,  which  is  assumed  stable  (often,  M  is 
finite- dimensional  LTI),  and  A  is  a  perturbation,  which  is  assumed  to  have  spatial  structure 

A  =  diag[8ilri,. .  .,6lITl,Al+i,- ■  -,Al+f\  (1) 

In  (1)  and  throughout  this  paper,  the  notation  diag[D\, . . . ,  Dn]  refers  to  a  block  diagonal 
matrix  with  blocks  Dk .  The  blocks  in  A  can  in  general  represent  real  parameters  or  dynamic  (LTI, 
LTV,  nonlinear)  perturbations. 

In  this  paper  we  consider  linear,  discrete-time  systems,  but  the  results  extend  with  minor 
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Figure  1:  Robust  stability 

changes  to  the  continuous  time  case.  The  I2  norm  is  used  for  signals:  '1S  the  set  of  Cn  valued 

square-summable  sequences  over  the  positive  integers.  Cc{f 2)  denotes  the  set  of  linear,  bounded 
and  causal  operators  in  Ilf.  All  our  perturbation  structures  will  therefore  be  a  subset  of  the  class  of 
structured  LTV  operators  of  the  form  (1)  with  6{  G  &L+j  G  T-dl™3),  1  <  i  <  L,  1  <  j  <  F. 

Some  of  the  blocks  will  additionally  be  specialized  to  be  LTI  operators,  or  further  to  be  real 
parameters.  The  uncertainty  is  normalized  to  a  ball  Ba  =  {A  :  ||A||  <  1}  in  the  /2-induced  norm. 

Assuming  M,  A  are  causal,  the  standard  notion  of  stability  of  the  interconnection  of  Figure 
1  is  that  the  map  between  injected  disturbances  d\,  d 2  and  the  signals  z,  w  is  causal,  and  it  is  a 
bounded  operator  when  restricted  to  I2.  If  M.  A  are  themselves  stable,  this  reduces  to  testing  for 
the  invertibility  of  /  —  AM.  This  is  captured  in  the  following: 

Definition  1  Assume  that  M  G  £(^2)-  The  system  of  Figure  1  is  robustly  stable  if  I- AM  :  /2 — ^^2 
has  a  causal,  bounded  inverse  for  all  A  E  Ba-  Robust  stability  is  uniform  if 

sup  || (7  —  AM)-1 1|  <  00  (2) 

AeBA 

This  paper  is  concerned  with  tests  for  robust  stability  under  various  assumptions  on  the  un¬ 
certain  perturbations.  An  important  comment  is  that  these  tests  can  also  be  used  for  analyzing 
robust  disturbance  rejection:  a  performance  specification  in  terms  of  a  bound  on  the  induced  norm 
of  a  transfer  function  can  under  fairly  general  circumstances  (see  [8,  18])  be  converted  to  a  robust 
stability  problem  with  an  additional  uncertainty  block. 
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2.1  Constant  Matrix  Analysis 

The  effect  of  the  structure  of  the  perturbation  in  the  robustness  analysis  is  apparent  when  we 
consider  the  constant  matrix  version  of  the  interconnection  of  Figure  1.  Now  M  ,  A  are  matrices 
in  CnX”,  A  still  has  the  spatial  structure  (1).  The  invertibility  of  I  —  AM  is  captured  by  the 
structured  singular  value  p  [5,  11,  7,  19],  defined  as  follows  (where  a  denotes  maximum  singular 
value): 

Definition  2  The  structured  singular  value  p  a  ( M )  of  a  matrix  M  with  respect  to  a  structure  A 
is  defined  as  p/\(M)  :=  0  if  no  A  £  A  makes  I  -  AM  singular,  otherwise 

(J-a(M)  :=  ( min  {d  (A)  :  det[7  —  AM]  =  0})-1  (3) 

Thus  /  -  AM  is  invertible  for  all  A  E  Ba  if  and  only  if  Pa(M)  <  1.  Two  different  cases  are 
of  interest:  complex  p  [11]  for  complex  A,  and  mixed  p  [19],  where  some  of  the  blocks  in  A  are 
restricted  to  be  real  (see  (8)  below).  In  this  paper  we  shall  use  a  common  generic  notation  p;  the 
distinction  will  be  made  explicit  whenever  is  necessary.  The  following  basic  property  of  p  is  known 
as  the  main  loop  theorem  [11]: 


Lemma  1 

as  M  :=  \ 


Given  a  block  structure  A  =  diag[ Ai,  A2]  and  a  C-valued  matrix,  suitably  partitioned 

M,U  1  then  pa(M)  <  1  if  and  only  if 
M21  M22 


Pa 2(^22)  <T  max  pAl(^*A2)<l 

a(A2)<i 

where  M  *  A2  is  the  Linear  Fractional  Transformation  (LFT) 


M  *  A2  :=  Mu  +  -Mx2A2(/  —  M22A2)  1M2i 


(4) 


Since  exact  computation  of  the  structured  singular  value  is  hard,  it  is  usually  approximated  by 
upper  and  lower  bounds.  For  lower  bound  algorithms  refer  to  [11,  19].  A  computationally  tractable 
upper  bound  is  obtained  by  considering  scaling  matrices  which  commute  with  the  elements  in  A. 
The  matrices  with  that  property  are  of  the  form 
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(5) 


X  —  diag  [Xi, ...  ?  X^,,  x ■■■■•,  x jJ^.pIrnF] , 

Let  X  be  the  set  of  positive  matrices  of  the  form  (5),  define 

A a(M):=  ud&iXMX-1)  (6) 

Then  pa (M)  <  pa(M).  Equivalently,  the  Linear  Matrix  Inequality  (LMI)  condition 

3X  E  X  :  M*XM  -  X  <  0  (7) 

is  sufficient  for  pa(M)  <  1.  This  convex  feasibility  condition  is  attractive  for  computation  (see 
[3]).  It  can  also  be  refined  for  the  case  of  mixed  uncertainty  structures.  If  A  is  of  the  form 

A  —  diag  \S\ITl , . . . ,  6lrI,  Sp^il . . . ,  ^Lr+Lc ?  ^L+h  •  •  •  >  (8) 

where  Si . .  .6lr  6  R  (i.e.  the  first  Lp  scalar  times  identity  blocks  correspond  to  real  parameters), 
then  condition  (7)  can  be  tightened  (see  [7,  19])  to 

3X  e  X,  G  6  Q  :  IkPXM  -  X  +  j(M*G  -  G*M)  <  0  (9) 

where  the  matrices  G  are  of  the  form  G  =  diag[G\, . . . ,  Gir,  0, . . . ,  0]  with  G;  =  G*. 

2.2  Robust  Stability  Tests 

With  the  notation  developed  from  the  constant  matrix  case,  we  now  summarize  known  conditions 
for  analysis  of  the  robust  stability  question  in  Figure  1,  where  M{erjJ)  is  assumed  to  be  always  a 
finite  dimensional  LTI  system. 

The  first  result  (see  [11])  is  that  if  A  is  LTI  uncertainty,  robust  stability  is  equivalent  to  the 
complex  p-test  across  frequency 

m&xpA{M{e^))  <  1  (10) 

LO 

An  analogous  condition  holds  if  the  uncertainty  structure  consists  of  a  combination  of  real  para¬ 
metric  and  LTI  perturbations,  the  only  difference  being  that  complex  p  is  replaced  by  mixed  p. 
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Convex  conditions  for  robust  stability  analysis  follow  from  the  //-upper  bound.  The  scalings  X 
can  be  chosen  either  to  be  constant,  or  frequency  varying,  giving  the  following  two  tests: 

inf||XM(e>)X-1||co<l  (11) 

ma x/ti\(M(e-7a;))  =  ^inf  j|X(w)M(eJ'a'l)X(u;)-1||oo  <  1  (12) 

Clearly,  (11)  implies  (12)  which  in  turn  implies  (10),  so  both  conditions  (11-12)  are  sufficient  for  the 
case  of  LTI  uncertainty.  Also,  for  mixed  real  parametric/LTI  perturbations,  (12)  can  be  tightened 
by  use  of  frequency  dependent  (S'-scales  as  in  (9). 

Recent  results  have  provided  an  exact  characterization  of  conditions  (11)  and  (12).  It  was  shown 
independently  by  Shamma  [18]  and  Megretski  [9]  that  the  constant  scales  test  (11)  is  necessary  and 
sufficient  for  robust  stability  under  structured  LTV  perturbations.  In  the  frequency  varying  case, 
Poolla  and  Tikku  [14]  have  recently  shown  that  (12)  is  necessary  and  sufficient  for  robust  stability 
against  the  class  of  arbitrarily  slowly  varying  structured  perturbations.  The  following  definition  is 
from  [14]. 

Definition  3  An  operator  A  £  £c(h)  has  time  variation  slower  than  u  if  ||AA—  AA||  <  z/||A||, 
where  A  is  the  delay  operator.  The  set  of  such  operators  is  denoted  by  T(v). 

The  norm  ||  A  A  -  AA||  is  a  natural  way  to  capture  the  rate  of  time  variation  of  an  operator.  If  this 
norm  is  zero,  A  commutes  with  A  and  the  operator  is  time  invariant.  Therefore  small  values  of  v 
in  Definition  3  correspond  to  operators  which  vary  “slowly”.  A  value  v  =  2  includes  all  “arbitrarily 
fast”  operators  on  1 2.  The  result  from  [14]  is  that  the  (12)  holds  if  and  only  if  the  system  is  robustly 
stable  under  perturbations  in  tF(o)  for  some  v  >  0. 

It  may  be  argued  from  these  results  that  //-analysis  should  be  abandoned  in  favor  of  these  upper 
bounds  which  appear  more  tractable,  especially  in  the  case  of  (12)  which  has  mild  conservatism. 
There  are  still  good  reasons,  however,  to  formulate  a  problem  in  terms  of  //.  In  the  first  place, 
although  the  upper  bounds  have  guaranteed  polynomial- time  computation,  the  size  of  the  problems 
can  be  very  large  and  render  the  computation  impractical.  In  these  cases  one  often  relies  heavily  on 
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the  availability  of  efficient  lower  bound  algorithms  [11,  19,  1]  (which  have  no  guarantees  but  appear 
to  behave  well  in  practice)  to  compute  the  analysis.  Secondly,  if  there  is  parametric  uncertainty 
in  the  problem,  the  upper  bounds  may  be  substantially  conservative  (there  is  no  corresponding 
slowly- varying  interpretation).  Lower  bound  algorithms  provide  a  fast  method  to  obtain  “bad” 
parameter  values,  and  can  be  further  employed  to  assess  this  conservatism  and,  if  desired,  pursue 
a  more  refined  analysis  by  branch  and  bound  techniques  [10]. 

2.3  A  Power  Distribution  Lemma 

The  following  Lemma  from  Poolla  and  Tikku  [14]  provides  a  useful  characterization  of  time  varying 
perturbations,  which  will  be  used  in  this  paper. 

Lemma  2  Let  0  <  oq  <  . . .  <  ur  <  tt  be  distinct  frequencies.  If  the  vector  valued  signals 

r  r 

z  =  Zke^Wkt,  w  —  Wk€^Ukt  (13) 

k=l  k— 1 

satisfy  the  power  inequality  J2k=i  ll^ll2  -  SL=i  INfc||2>  then  there  exists  a  linear  time-varying, 
causal  operator  A  such  that 

(i)  || A||  <  1 

(ii)  A  £  F(v),  v  =  2sin(a,r  i°l) 

(Hi)  A z  =  w  +  wtr ,  wtr  G  h 


Heuristically,  this  lemma  says  that  provided  that  the  total  power  of  z  is  greater  than  that  of  w,  a 
contractive  LTV  operator  can  rearrange  the  power  between  frequencies,  mapping  z  to  w  in  steady 
state.  The  time  variation  v  required  is  a  function  of  the  amount  of  “frequency  shifting”  performed. 
In  contrast,  a  contractive  LTI  operator  will  always  decrease  the  power  at  every  frequency. 
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3  A  //-Test  for  Mixed  LTV/LTI  Analysis 

Sections  3  and  4  refer  to  the  system  of  Figure  1,  where  we  now  set 


Mg $  Mg  _  ’ 


(14) 


In  (14),  M  is  a  finite  dimensional  LTI  discrete  time  system,  which  can  also  be  given  a  state-space 
representation 

Equivalently,  M  is  obtained  from  the  LFT  M( A)  =  XI  *  S  =  D  +  C\(I  —  XA)~1B,  where  A  is  the 
delay  operator  (or  a  frequency  variable  in  the  unit  disk)  and  S  is  the  constant  state-space  matrix. 

$  is  a  causal  LTV  operator  on  1™,  and  0  is  a  causal  LTI  operator  on  l9-  Each  has  a  spatial 
structure  analogous  to  (1),  denoted  respectively  by  4>  and  0.  B<$  and  B©  denote  the  unit  balls  of 
uncertainty.  X4,  X9  will  denote  the  sets  of  scaling  matrices  corresponding  to  each  structure.  The 
consideration  of  real  parametric  uncertainty  is  deferred  to  Section  5. 

An  important  integer  parameter  determined  by  the  structure  is  the  dimension  d  of  the  space  of 
hermitian  scaling  matrices  which  commute  with  the  LTV  structure  4?.  (X®  is  the  positive  cone  in 
this  space).  When  consists  only  of  full  blocks,  d  is  the  number  of  blocks. 

The  main  result  in  this  section  is  an  extension  of  the  exact  //-test  (10),  for  LTI  analysis,  to  the 
case  of  mixed  LTV/LTI  perturbations. 

Since  LTV  structures  are  usually  characterized  by  tests  in  terms  of  ft,  rather  than  ft,  it  is  not 
obvious  that  a  ft- test  can  capture  the  mixed  LTV/LTI  case  (and  in  particular,  the  LTV  case).  The 
main  idea  to  obtain  this  //-test  is  inspired  in  work  by  Bercovici  et.al.[2],  where  an  augmentation  or 
lifting  in  the  structure  converts  the  upper  bound  ft  to  ft  of  a  larger  matrix.  The  results  in  [2]  apply 
to  constant  matrices  and  are  based  on  operator-theoretic  methods,  but  can  also  be  obtained  as  a 
corollary  of  the  more  general  dynamic  results  to  be  presented  in  this  paper,  which  will  be  proved 
by  convex  analysis  methods. 
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We  first  consider  the  augmented  matrix 


M(wi,  ...,ud)  =  diag{M(e3W1), . . .,  M{e*u*)\ 


(16) 


which  amounts  to  considering  the  system  M  at  a  fixed  number  d  of  frequencies.  Next  we  introduce 
the  following  augmented  structures  in  the  space  of  complex  matrices: 


$  = 


A  = 


- 1 

(-1 

fydd. 

'  A11  . 

.  .  Ald  ' 

Adl  . 

/\dd 

, G  i,j  =  1  ...d 

(  A”  G  A,  i  =  1 . .  .d 
’  1  0],  i  ±  j 


(17) 


(18) 


The  augmented  structures  <£,  A,  are  d  times  larger  than  the  corresponding  d>.  A.  For  the  case 
of  <&,  it  is  obtained  simply  by  considering  d2  copies  of  #,  in  matrix  form.  A,  which  contains  $ 
as  a  submatrix,  is  obtained  in  a  similar  fashion,  the  only  difference  being  that  the  time  invariant 
blocks  6  are  only  “copied”  along  the  diagonal,  and  the  rest  of  the  entries  are  set  to  zero.  As  an 
illustration,  Figure  2  (a)  contains  the  augmented  configuration  (M,  A)  for  the  case  d  =  2. 


(«)  (b) 


Figure  2:  Augmented  representations 
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This  structure  has  a  “frequency  shifting”  interpretation  which  relates  to  the  remarks  made 
in  regard  to  Lemma  2:  the  augmentation  corresponds  to  considering  a  system  at  a  number  d  of 
frequencies,  and  the  different  treatment  of  LTI  and  LTV  blocks  is  due  to  the  fact  that  only  the 
time- varying  perturbations  are  allowed  to  “shift  energy  between  frequencies”;  this  is  represented  in 
$  by  the  off-diagonal  terms. 

It  is  also  convenient  to  consider  the  configuration  of  Figure  2  (6),  where  the  LTI  portion  is 
included  with  the  dynamics.  Define 

G(u,  6)  :=  M{ejw)  *6  =  M$  +  M<s>e0(I  -  Me0y1Me<s,  (19) 


For  this  LFT  to  be  well  defined  for  6  6  B©,  the  following  condition  must  hold: 

ma x^©(M^(e7U'))  <  1  (20) 

(jJ 

Under  this  condition,  G(ui,6)  is  continuous  for  u>  G  [— 7r,7r],  9  G  B©.  For  given  (ajk,6k),k  =  1 . .  .d, 
the  matrix  G(u>i, . . .  ,u>^,  91 , . . .,  6d)  :=  diag[G(u> i,  61) . . . G(u>d ,  6d)\  is  obtained  by  LFT  between  the 
9  portion  of  A  and  the  matrix  . . .  ,u>d ),  as  shown  in  Figure  2  ( b ). 

We  now  state  the  main  result: 


Theorem  3  In  reference  to  the  system  (14),  the  following  are  equivalent: 

(a)  The  system  is  uniformly  robustly  stable. 

(b)  Condition  (20)  holds,  and  with  u>k  varying  in  [-7r,7r],  and  6k  G  B©, 

max  ni(G(oJi,...,u)d,91,...,ed))  <  1  (21) 

(c)  With  u>k  varying  in  [— 7r,  7t] 


max  <  1 

LUl,...,Ud 


(22) 


Proof:  The  equivalence  of  (6)  and  (c)  is  a  simple  consequence  of  the  main  loop  theorem  (Lemma 
1).  In  fact,  for  each  wi, . . .  Lemma  1  implies  that 


|  /XQ(Me(ejwk))  <  1,  k  =  1 . .  .d 

l  max<?i  -,&*))  <  1 


(23) 
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Taking  maximum  over  wi, . . .  ,u>d,  the  right  hand  side  of  (23)  gives  (20)  and  (21). 

The  equivalence  of  (a)  and  (6)  follows  from  a  more  technical  convex  analysis  argument,  and  is 
covered  in  the  Appendix. 

Theorem  3  provides  a  necessary  and  sufficient  test  for  robust  stability  in  terms  of  a  //-condition 
(22)  which  involves  a  search  in  d  frequency  variables.  Using  a  state-space  realization  as  in  (15),  a 
state-space  test  in  terms  of  g  of  a  constant  matrix  can  be  derived,  in  an  analogous  way  as  in  the 
standard  LTI  problems  involving  one  frequency  variable  ([11]): 


Corollary  4  The  system  has  uniform  robust  stability  if  and  only  if  <  1,  where  S  = 

diag[S, . . .,  S'],  and 


As  = 


(24) 


A  I1  ...  Ajf 


A|  =  diag[XI,  6%  £  $,  i=l...d 

A5  =  diag[0,  4>*J,  0],  §13  £  i  ^  j,  i,j  =  l...d 


(25) 


Proof:  M{ A)  (with  A  a  complex  variable  in  the  closed  unit  disk  D),  is  the  transfer  function  of 
the  stable  state-space  system  given  by  S.  The  main  loop  theorem  (Lemma  1)  implies  that 

/^As(^)  <  1  ^  max  _fiA(M(Xu...,Xd))  <  1  (26) 

Ai,...,A  d£D 

where  M(X\, . . . ,  Xd)  '•=  diag[M(Xf), . .  .,M(A</)].  The  maximum  modulus-like  property  [11]  of 
complex  g  implies  that  the  maximum  in  (26)  occurs  at  the  boundary  of  the  disk,  and  therefore  (26) 
coincides  with  (c)  in  Theorem  3. 


□ 
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4  Convex  Tests  for  the  Mixed  LTV /LTI  Problem 

The  previous  tests  have  the  theoretical  advantage  of  being  exact,  but  as  remarked  before,  exact  [i 
computation  is  hard,  so  practical  use  of  these  conditions  will  involve  employing  bounds  as  those  in 
standard  software  packages  such  as  /x- Tools  [1]. 

In  particular,  we  analyze  in  this  section  the  upper  bounds  for  this  problem  (sufficient  conditions 
for  robust  stability)  which  lead  to  convex  optimization.  For  this  purpose,  the  g  conditions  obtained 
in  Section  3  can  be  bounded  by  use  of  the  constant  matrix  upper  bound  /x.  We  focus  here  on  the 
upper  bound  over  the  frequencies  uq,. . .  which  follows  from  (22).  An  examination  of 

the  augmented  structure  A  shows  that  the  corresponding  commuting  matrices  X  are  of  the  form 

1  =  diag[X*, Xl  Xl, . . . ,  Xed]  (27) 

where  6  X$,  X{  6  X*. 

An  alternative  is  to  directly  apply  scaled  small-gain  conditions  to  the  original  problem.  In  the 
case  of  mixed  LTV/LTI  analysis  as  in  (14),  the  natural  scaling  set  is  of  the  “mixed”  form 

X(oj)  =  diag[X $,  Xe{u)\  (28) 

where  the  portion  X6  which  corresponds  to  the  LTI  blocks  6  is  allowed  to  vary  in  frequency,  and 
the  portion  X®  corresponding  to  the  LTV  blocks  is  constant.  Xm  will  denote  the  set  of  such  scaling 
functions;  without  loss  of  generality  they  are  assumed  to  be  continuous  over  frequency. 

We  now  show  that  the  two  approaches  are  equivalent. 

Proposition  5  Given  the  system  (14),  the  following  are  equivalent: 

max  .  ..,Wd))  <  1  (29) 

3  X(u)  6  xm  :  o  (W(u;)M(eiw)X(u;)-1)  <  1  Vu;  (30) 
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Proof:  If  an  X  satisfying  (30)  is  found,  then  for  any  choice  of  frequencies  uq, . . . ,  u>d,  setting 


1  =  diag[ X*t  X  V),  •  •  • ,  X\oJd)} 

will  result  in  a  (XM(uq, . .  .,u^)X_1)  <  1,  and  therefore  <  1,  implying  (29). 

The  converse  implication  is  covered  in  the  Appendix. 

As  a  corollary  of  Proposition  5  and  Theorem  3,  condition  (30)  is  sufficient  for  robust  stability 
under  mixed  LTV/LTI  perturbations;  this  could  also  be  shown  via  standard  small-gain  arguments. 

From  a  computational  point  of  view,  (30)  has  the  advantage  of  involving  a  search  over  only  one 
frequency  variable.  A  direct  approach  1  would  be  to  grid  the  frequency  axis  and  convert  (30)  to  an 
LMI  condition.  Note,  however,  that  the  common  scale  introduces  a  coupling  in  the  problem,  so 
one  is  left  with  a  large  LMI  condition,  with  size  growing  with  the  number  of  frequency  grid-points. 

In  comparison,  (29)  tells  us  that  in  fact,  d  frequency  values  suffice,  and  we  must  only  solve  a 
coupled  LMI  problem  of  this  size.  However,  since  these  frequencies  are  not  known  a  priori,  one  has 
to  grid  a  d  dimensional  space  of  frequencies.  Therefore  (29)  reduces  the  size  of  the  coupled  LMI  at 
the  expense  of  more  gridding.  For  low  values  of  d.  this  alternative  may  be  convenient. 

We  now  consider  the  question  of  the  conservatism  of  these  conditions.  The  results  of  [14] 
reviewed  in  Section  2  suggest  that  the  conditions  become  necessary  if  the  LTI  perturbations  are 
enlarged  to  include  arbitrarily  slowly  varying  uncertainty.  This  is  the  content  of  the  following 
statement;  the  proof  involves  an  extension  of  the  techniques  of  [14],  described  in  the  Appendix. 

Theorem  6  The  conditions  in  Proposition  5  are  satisfied  if  and  only  if  there  exists  u  >  0  such  that 
the  system  (14)  has  uniform  robust  stability  for  A  =  diag[§,6],  $  arbitrary  structured  operator,  6 
structured  operator  in  X(v). 

1  Another  approach  would  be  to  parametrize  X6{io)  by  basis  functions. 
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5  Combination  with  Real  Parametric  Uncertainty 


In  this  section  we  take  a  further  step  in  the  analysis  under  combined  uncertainty  structures;  in 
addition  to  LTV  and  LTI  blocks,  we  include  real  parametric  perturbations.  We  will  consider  the 


robustness  analysis  setup  of  Figure  3: 


Figure  3:  Robust  stability  problem 


$  = 


$00 
0  0  0 

0  0  <5 


A  0  ' 
0  6 


H  = 


//$  H<t>o 
H$<t>  He  Hgs 
Use  Hs 


'  H a  Has  ' 
.  Hsa  Hs 


(31) 


In  (31),  H{eiw)  is  a  stable  system.  $  and  0  are  structured  LTV  and  LTI  perturbations  as  before. 
A  is  defined  as  in  (14).  The  additional  structured  perturbation  S  consists  of  real  parametric  blocks 
(e.g.  6  =  diag[S\I , . .  .,SnI],  Si  €  R).  The  notation  ($,  H  instead  of  A,  M)  is  chosen  to  clarify  the 
proofs  below. 

We  wish  to  obtain  a  necessary  and  sufficient  condition  for  robust  stability  in  this  class,  extending 
the  results  in  Section  3.  Since  real  parameters  are  a  special  case  of  time  invariant  uncertainty, 
at  first  sight  it  would  appear  that  Theorem  3  applies  directly,  yielding  a  m i xed-/i  condition  on 
H  —  diag[H(e3Wl), . . . ,  H(e?Wd)],  analogous  to  (22)  but  with  the  copies  of  S  in  the  augmented 
structure  constrained  to  be  real. 

This  augmentation  would  not  capture,  however,  an  additional  property  of  real  parametric  uncer¬ 
tainty:  in  addition  to  taking  real  (rather  than  complex)  values,  a  real  parameter  has  no  dynamics 
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and  therefore  is  constrained  to  be  constant  across  frequency2 .  This  suggests  a  modification  of 
the  augmented  uncertainty  structure  for  the  case  of  real  parameters,  where  they  are  forced  to  be 
constant  across  the  augmentation.  Consider  the  structure 


§  = 


tp11 
i Pdl 


...  ' 
. . .  1pdd  _ 


•{ 


ipu  =  diag[$n ,  9l,6],i  =  1 . .  .d 
=  diag[$lj ,  0,  0],  i  ^  j 


(32) 


where  £  <I>  and  9l  £  0  as  before,  but  we  constrain  the  copies  of  6  to  be  repeated  across  the 
augmentation.  The  d  =  2  case  is  depicted  in  Figure  4  (a). 


(a)  ( b ) 

Figure  4:  Augmented  systems  for  the  real  parametric  case 


It  will  also  be  useful  to  “close  the  loop”  on  the  real  parametric  part,  which  will  reduce  the 

2This  fact  does  not  come  into  play  in  standard  mixed  /r  with  LTI/parametric  uncertainty,  where  only  one  frequency 
is  involved  in  the  destabilization. 
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problem  to  the  situation  of  Section  3.  Assume  that  the  real  p,  condition 


ma <  1  (33) 

(jJ 

holds,  then  for  any  fixed  6  6  B<$  we  can  define 

M(ejw,6)  :=  H{C“)  *6  =  HA  +  HAS6{I  -  HS6)~1HSA  (34) 


which  is  a  stable  time  invariant  system  for  each  6  6  B§. 

For  fixed  6,  and  given  uq  ,...,Ud,  the  matrix  M(uq, . .  .,ud,S)  :=  diag[M{e^1,b) . . .  M(e^d,  6)] 
(which  corresponds  exactly  to  the  augmentation  (16))  is  obtained  by  LFT  between  the  6  portion 
of  IP  and  the  matrix  H(oj  1, . .  ,,c jj).  The  uncertainty  structure  corresponding  to  M  is  A  as  in  (18), 
which  is  depicted  in  Figure  4  ( b )  for  the  case  d  =  2. 

We  now  provide  the  extension  of  Theorem  3. 


Theorem  7  In  reference  to  the  system  (31),  the  following  are  equivalent: 

(a)  The  system  is  uniformly  robustly  stable. 

(b)  Condition  (S3)  holds ,  and  with  u>k  varying  in  [— 7r,7r],  6  £  B5 

max  ..,Ud,6))  <  1  (35) 

wi, ...,LOd,  S 

(c)  With  ujk  varying  in  [— 7T,7r], 

max  ns.(H{uu...,ud))  <  1  (36) 

Ull,...,Uld 

where  /i  is  the  mixed  (real/ complex)  structured  singular  value  with  respect  to  the  structure  (32). 


Proof:  The  equivalence  of  ( b )  and  (c)  is  a  direct  application  of  Lemma  1. 

Also,  (a)  implies  condition  (33),  otherwise  the  real  parameters  S  would  destabilize  by  themselves. 
Given  (33),  we  observe  that  for  fixed  6  £  Bj, 


I-  qH  = 


[  /  -AHas(I  -  SHs)-1  ' 

'  /  -  A(H  *  S)  0 

[  0  I 

-SHs A  I  -  6H5 

(37) 
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It  follows  that  the  (uniform)  robust  stability  of  system  (31)  is  equivalent  to  the  fact  that  for 
every  b  in  B5,  the  system  6),  A)  is  (uniformly)  robustly  stable.  Since  A  has  the  mixed 

LTV/LTI  structure  of  Theorem  3,  this  is  in  turn  equivalent  to  condition  (35)  in  (6). 

□ 

The  previous  result  has  reduced  the  robust  stability  problem  under  LTV,  LTI  and  parametric 
uncertainties  to  a  mixed-/i  condition  across  d-frequencies.  We  remark  the  following: 

•  In  analogous  manner  to  Corollary  4,  a  state-space  condition  can  be  derived  from  (36),  which 
is  equivalent  to  a  single  mixed  fi  problem. 

•  As  usual,  practical  computation  of  a  ^-condition  such  as  (36)  must  be  approached  by  upper 
and  lower  bounds,  and  possibly  branch  and  bound  techniques.  Upper  bounds  will  have  the 
form  given  in  (9);  analogously  to  the  situation  in  Section  4,  we  have  the  choice  of  writing 
these  conditions  in  the  original  problem  or  in  the  augmented  problem.  In  this  case,  however, 
these  are  not  equivalent:  since  the  augmentation  introduces  repetition  in  the  real  uncertainty, 
this  increases  the  freedom  of  the  scaling  matrices  Xs ,  Gs  corresponding  to  these  blocks, 
which  provides  a  way  of  imposing  the  condition  that  S  is  constant  across  frequency.  Without 
augmentation,  this  condition  is  not  imposed  so  the  bound  is  weaker. 


6  Examples 

In  this  section  we  illustrate  the  results  of  this  paper  with  a  series  of  examples.  These  have  been 
deliberately  constructed  so  that  there  is  a  direct  way  to  answer  and  interpret  the  robust  stability 
question,  thus  providing  more  insight  into  the  conditions  given  in  the  previous  sections. 

6.1  A  system  with  LTV  uncertainty 

We  consider  the  interconnection  of  Figure  (5),  where  Gi,  G2  are  LTI  SISO  systems,  and  $1,  $2 
are  uncertain  perturbations  with  ||4q||  <  1. 
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Figure  5:  Example  of  analysis  of  LTV  uncertainty 


•  If  $1,  $2  were  LTI  perturbations,  then  they  would  commute  with  G\,  G 2  and  from  the 
small-gain  theorem  the  robust  stability  condition  would  be 


||<SfiGr2||00  <  1 


(38) 


•  From  now  on  we  consider  the  more  interesting  case  when  $2  are  arbitrary  LTV  pertur¬ 
bations.  It  turns  out  that  in  this  simple  configuration,  the  necessary  and  sufficient  condition 
for  robust  stability  is 

I|GiIU|G2|U  <  1  (39) 

This  condition  is  in  general  stricter  than  (38),  since  the  two  transfer  functions  G\{e^w)  and 
need  not  achieve  their  peak  gain  at  the  same  frequency. 


The  sufficiency  of  (39)  is  clear  from  small-gain.  To  explain  why  it  is  also  necessary,  let  us 
particularize  in  the  example  (where  A  is  the  delay) 


Gr 


0T 

1  -  0.9A’ 


G-2 


-0.1 
1  +  0.9A 


(40) 


Here  both  systems  have  H0 o  norm  equal  to  one,  achieved  respectively  at  frequencies  0, 7 r.  Their 
magnitude  frequency  response  plots  are  depicted  in  Figure  6  below.  Condition  (39)  is  therefore  not 
satisfied.  Let  us  show  that  the  system  can  be  destabilized  by  contractive  LTV  operators.  In  this 
case  it  suffices  to  consider  the  time-varying  gains  3>i(f)  =  $2 (t)  =  (— 1)<- 
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In  reference  to  Figure  5,  consider  a  constant  signal  v(t)  =  1  at  the  input  of  G\.  Since  the  transfer 
function  G\  is  1  at  u>  =  0,  the  steady  state  output  will  be  the  same  signal.  The  time-varying  gain 
$i(t)  modulates  this  signal  to  w(t)  =  (-1)*,  which  has  all  its  frequency  content  at  to  =  7r,  where 
G2  has  value  1.  This  implies  the  steady  state  output  of  G2  is  w(t),  which  is  demodulated  back 
to  v(t)  by  $2 (t).  We  have  a  steady  state  signal  in  the  loop,  which  violates  robust  stability.  This 
informal  argument  can  be  formalized  and  extended  to  arbitrary  G i,  G 2,  and  it  illustrates  strongly 
the  “frequency-shifting”  properties  of  LTV  perturbations. 


We  should  recover  the  same  answer  if  we  do  the  analysis  using  the  results  in  Section  3.  For  this 
purpose,  we  first  rearrange  Figure  5  to  an  M- A  setup,  where  A  contains  in  this  case  only  the  LTV 


portion  $ 


A  =  $  = 


>1  0  ' 

0  $2 


0  G\ 
G2  0 


Since  d  =  2,  we  must  compute  the  augmentation  of  Theorem  3  over  two  frequency  variables  uq, 
u2.  Figure  7  contains  the  resulting  function  n^(M(u>  1,^2)),  computed  using  the  software  package 
/u-Tools  [1].  We  find  that  the  maximum  is  1,  achieved  when  the  two  frequencies  take  the  values  0,  w, 
which  is  consistent  with  the  previous  analysis.  Similar  results  can  be  obtained  using  the  state-space 
condition  in  Corollary  4. 
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Figure  7:  as  a  function  of  wi,  u2 

6.2  Including  LTI  uncertainty 

We  now  modify  the  previous  setup  to  include  LTI  uncertainty;  Gi,  G2  are  replaced  by 

P  _ .  r<  _ _ .  im.ii  <  i 

Gl  1  —  (0.6  +  0.3#i)A’  2  1  + (0.6- 0.3  W  11  4,1  “ 

which  are  LFTs  Gi  =  T\  *  0,  on  the  LTI  perturbations  6n  as  depicted  in  Figure  8. 


Figure  8:  System  with  LTV  and  LTI  uncertainty 


(41) 
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From  (39),  it  is  clear  that  the  worst-case  0,-  are  those  which  achieve  max^  The  values  in 

(41)  have  been  chosen  so  that  the  worst-case  perturbations  are  Q\  =  l,  02  =  —1  ,  which  produce  G; 
as  in  (40).  From  the  previous  analysis  the  smallest  destabilizing  perturbation  is  of  norm  1:  6 j  =  1, 
62  =  —1,  $1  =  $2  =  (— 1)*. 

These  results  are  verified  when  we  do  the  analysis  of  Theorem  3.  After  rearranging  Figure  8 
in  the  standard  setup,  and  performing  the  augmentation,  we  obtain  max^ <W2  fiA(M(oJi,u>2))  =  1 
(achieved  at  uq  =  0,  u2  =  x). 

If  the  6i  in  (41)  correspond  to  real  parameters  instead  of  LTI  perturbations,  we  still  obtain  the 
same  answer  from  the  robust  stability  analysis,  since  the  worst-case  values  obtained  above  happen 
to  be  real.  This  can  also  be  verified  by  computing  the  mixed- /i  condition  on  the  augmented  system. 


6.3  Real  parametric  vs  LTI  perturbations 

To  produce  an  example  where  real  and  LTI  perturbations  give  a  different  answer,  one  can  use  the 
same  structure  as  in  (41),  but  impose  #1  =  02  =  9  (repeated  perturbation).  In  other  words, 


1  1  —  (0.6  +  O.30)A’  2  1  +  (0.6  -  O.30)A  v  ; 

If  6  is  LTI,  the  repetition  does  not  alter  the  results  in  Section  6.2,  since  an  LTI  perturbation  can 
take  the  values  1  at  to  =  0  and  — 1  at  m  =  ir.  An  example  is  9(eJUJ)  —  e3UJ  ( 6  is  the  delay  operator 
A),  which  turns  the  Gi  into  second  order  systems  with  IIGiH^  =  Halloo  =  1. 

The  same  answer  is  obtained  from  the  augmented  /./,  test  (22)  which  gives 

max  fix  (M(u>i,lo2))  =  1  (43) 

Ul-1 ,0/2 

achieved  at  =  0,  ui2  =  7T,  and  a  destabilizing  perturbation  with  01  =  1,  02  =  — 1,  as  expected. 

We  now  change  0  for  a  real  parameter  6  £  [-1, 1],  repeated  in  Gi,  G2 .  Since  it  is  constant  across 
frequency,  it  cannot  maximize  both  simultaneously.  In  fact,  straightforward  calculations 


show  that 


ll^llloo  Halloo  - 


16  -  9<52 
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(44) 


which  has  a  maximum  of  y  for  6  G  [—1, 1].  This  implies  by  (39)  that  the  overall  system  is  stable 
for  ||$, -||  <1,6  6  [-1, 11- 

Proceeding  by  augmentation,  (36)  gives  a  value 


maxui(j(wi,U2))  =  0.7906  <  1 

,U/2 


(45) 


In  this  case  we  can  show  directly  that  the  smallest  destabilizing  perturbation  for  6  6  R  and 
LTV,  has  norm  0  7g06  =  1.2649.  For  this  purpose,  choose  6  £  [—7,7]  and  ||$t-||  <  7.  The  full 
system  will  be  stable  as  long  as 


1 


16-97a=«n-r](l|G,IUI|GjlL)< 

Condition  (46)  is  equivalent  to  7  <  y/8/5  =  1.2649,  as  expected. 


(46) 


7  Conclusions 

This  paper  shows  that  a  combination  of  different  classes  of  uncertain  perturbations  (LTV,LTI,  para¬ 
metric)  can  be  analyzed  with  the  same  mathematical  tools  as  non-mixed  problems.  A  structured 
singular  value  condition  was  obtained,  applicable  to  any  combination  of  these  uncertainty  classes. 

These  results  allow  a  totally  decoupled  approach  to  uncertainty  modeling  in  complex  systems: 
one  can  choose  the  most  adequate  uncertainty  description  (LTV,  LTI,  parametric)  at  the  subsystem 
level,  and  obtain  an  exact  condition  for  robustness  analysis  of  the  overall  system  under  the  combined 
uncertainty  structure.  The  price  paid  in  terms  of  complexity  of  these  conditions  is  the  size  of  the 
corresponding  augmentation. 

From  a  computational  point  of  view,  a  number  of  equivalent  conditions  have  been  obtained, 
and  further  research  is  required  to  determine  the  most  efficient  approach  for  practical  problems. 
In  regard  to  the  convex  upper  bounds  which  lead  to  coupled  LMI  problems  across  frequency,  two 
alternatives  (29)  and  (30)  have  been  discussed  and  should  be  further  explored.  In  relation  to  lower 
bound  computation,  the  repeated  structure  of  the  augmented  systems  may  be  exploited  to  improve 
the  algorithms. 
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Finally,  scaled  small-gain  conditions  such  as  (30)  can  naturally  lead  to  the  extension  of  “D-K 
iteration”  methods  for  controller  synthesis.  In  the  author’s  opinion,  however,  synthesis  should  be 
based  on  simpler  and  more  heuristic  methods,  and  these  tools  with  very  specialized  uncertainty 
structures  are  best  employed  for  analysis  validating  the  resulting  designs. 


Appendix 

Convex  Analysis  Lemmas 

The  following  results  from  convex  analysis  will  be  used  in  the  proofs. 

For  a  set  K  C  Mf/,  co(K )  will  denote  the  convex  hull  of  K,  i.e.  the  set  of  all  convex  combinations 
of  elements  of  I{.  If  K  is  compact,  so  is  co(K). 

Lemma  8  (Helly)  Let  {KujweO.  be  family  of  convex  closed  sets  in  Rd,  of  which  at  least  one  is 
bounded.  If  Ku  =  0  then  there  exist  d  +  1  sets  KUJt ,  i  =  1 . . .  d  +  1  with  empty  intersection. 

Lemma  9  Let  K ,L  be  disjoint  convex  sets  in  RR  where  K  is  compact  and  L  is  closed.  Then  there 
exists  a  vector  x  £  Rd,  and  a,  /3  in  R  such  that 

(x,k)  <  a  <  /3  <  (x,l)  Mk  £  K,l  £  L.  (47) 

Lemma  10  If  K  C  Rd,  every  point  in  co(K )  is  a  convex  combination  of  d  +  1  points  in  K;  for  K 
compact,  every  point  in  the  boundary  of  co(K)  is  a  convex  combination  of  d  points  in  K . 

Proofs:  Lemmas  8  and  9  can  be  found  in  [16].  The  first  part  of  Lemma  10  is  a  result  of 
Caratheodory  (see  [16]),  and  implies  that  for  every  v  £  co(K),  there  exists  a  simplex  of  the  form 

(d+ 1  d+ 1  'l 

S(v i,...,vd+1)  =  j  Y^akvk  ■  otk  >  0,^0*  =  1  > 

lfc=l  k— 1  J 

with  vertices  Vk  6  K ,  which  contains  v.  If  the  Vk  are  in  a  lower  dimensional  hyperplane,  then 
d  points  will  suffice  to  generate  v.  If  not,  then  every  point  in  S(vx, . . . ,  ^-n)  corresponding  to 
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atk  >  0  VA:  will  be  interior  to  S(v i, . . . ,  vd+i)  C  co(K).  Therefore  for  points  v  in  the  boundary  of 
co(K),  one  of  the  a^s  must  be  0  and  a  convex  combination  of  d  points  will  suffice,  completing 
Lemma  10. 


Proof  of  Theorem  3 

We  now  show  the  equivalence  of  (a)  and  (6)  in  Theorem  3.  For  simplicity  we  will  write  the  proof 
for  the  case  where  the  structure  $  contains  only  full  blocks,  i.e.  $  =  diag[$> i, . .  .$d],  €  C(l™1). 

The  same  arguments  can  be  extended  to  the  case  of  61  blocks  in  the  uncertainty  in  a  similar 
style  as  in  [12].  The  following  Lemma  gives  an  interpretation  of  the  uncertainty  structure  $. 


Lemma  11  Assume  (20)  holds.  Given  uq, . .  ,,ud,  91,. .  .,0d  the  following  are  equivalent: 

(i)  1,...,wd,61,...,0d))  >  1. 

(ii)  3  v1,..  ,,vd  £  Cm,  not  all  zero  such  that 

£|| (G(uk,0k)vk),  2>^||(A  \  l  =  l...d 


fe=l 


fc= 1 


(48) 


Proof:  (?)  is  equivalent  to  the  existence  of  cr($)  <  1  and  a  vector  v  =  co/(u1, . .  .,vd)  ^  0 
such  that  (/  —  if>G(u i , . . . ,  u>d ,  91 , . . . ,  0d))v  =  0,  or 


v 1 

(2(aq,  91)v1 

vd 

. .  '■a 

.  G(u>d,6d)vd  _ 

Each  has  in  turn  d  subblocks,  which  will  henceforth  indexed  by  the  subindex  /,  and  impose 
a  partition  on  the  vk ,  G{ u>k,9k)vk.  The  equations  and  variables  in  (49)  can  be  reordered,  to  put 
together  the  vj, . . . ,  vf  and  (<7(aq,  01)u1)/, . . . ,  (G(ud,  6d)vd)i  for  each  /,  which  reduces  the  structure 
of  $  to  d  full  blocks  in  this  new  order.  Now  a(4>)  <  1  is  equivalent  to  the  norm  inequalities  (48). 


□ 
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(«)=►(&) 

Condition  (20)  is  clearly  necessary  for  robust  stability;  if  it  did  not  hold,  the  standard  results 
(10)  imply  that  the  system  could  be  destabilized  by  LTI  perturbations  alone.  Therefore  G{u,9)  is 
well  defined  and  continuous. 

Assume  there  exist  u\,. .  .,u>d  £  [— x,  7r],  01, . .  .,0d  E  B©  with  . . . ,  u>d,  91, . . . ,  9d ))  >  1. 

We  can  use  Lemma  11  to  obtain  v1, . . . ,  vd  satisfying  (48). 

Fix  e  >  0;  by  continuity  of  G  we  can  perturb  the  91 , . . . ,  9d  to  have  a(9k)  <  1  (strict  inequality), 
and  the  u\,. .  .,u>d  to  make  them  distinct,  satisfying 

£  ||(<?K,  9k)v%  ||2  >  (1  -  e)2  £  ||(  vk)t\\2  ,  l  =  l...d  (50) 

k= 1  k= 1 

We  are  now  in  the  conditions  of  an  interpolation  result  given  in  [6],  which  states  that  there  exists 

a  causal,  stable,  rational  LTI  perturbation  9{e^)  E  B©  3,  satisfying  9(e^Wk )  =  9k .  Now  define 

d  d 

v(t)  =  £  vke^ht,  z{t)  =  Y^(.G{uk,9k)vk)eiuJkt  (51) 

k—1  k= 1 

Then  G(u>,  9(e3U1)))v(t)  =  z(t)  +  e(t),  where  e(t)  is  a  transient  term.  For  any  /,  (50)  implies  that 
the  power  of  the  /-1th  component  of  z(t),  is  greater  than  the  power  of  the  /-1th  component  of  v(t), 
times  (1  —  e).  We  can  invoke  Lemma  2  to  show  the  existence  of  a  causal  time- varying  operator 
$1  which  maps  z(t)i  to  (1  —  e)v(t)i,  up  to  a  transient  term.  Constructing  $  =  diag[§\ . . .  we 
obtain 

(/  -  $G(w,  9(e3UJ)))v(t)  =  ev(t)  +  e(t)  (52) 


where  e(/)  E  h  is  transient.  This  implies  that 


(/  -  $(M  *  9))-1 1|  =  ||(7  -  $G)_1||  >  - 


Referring  to  (14),  it  can  be  verified  that 


I- AM 


'/  -  dMg)-1  r  7-  $(M*9)  0 

.  0  7  -9M9$  I  -  9Me 


3 [6]  constructs  a  perturbation  in  the  disk  algebra  »4(T)  interpolating  a  countable  number  of  frequencies;  in  the 
case  of  a  finite  number  of  frequencies  it  can  be  chosen  to  be  rational 
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This  implies  that 


||(/-  AM)-1!!  > 


f  I  -  $(Af  *6)  0 

-l 

—6Mq§  I  —  0M$ 

[  I  -  OMeY1 

0  / 

(55) 


As  e— >0,  the  numerator  in  (55)  tends  to  infinity  from  (53),  and  the  denominator  is  bounded  using 
(20).  This  implies  that  supA  ||(/  -  AM)-1||  =  oo,  violating  uniform  robust  stability. 


(*)=K«): 

For  this  direction  we  introduce  some  notation.  For  ui  £  [— 7r,  x],  0  £  Bq  C  CpXp,  v  £  Cm,  define 


<*(*)  =  ||(G(u;^)u)/||2  -  |H|2  ,  1=1...  d  (56) 

KAV)  =  (<£,«(»)>  •  •  •  >  atAv))  e  Kd  (57) 

V  =  {Aw,*(u)  :  w  £  [-x, x],  6  £  B0,  t>  £  Cm,  ||u||  =  1}  (58) 


From  (20),  AWte(v)  is  continuous  in  its  3  variables,  therefore  V  is  compact  in  Rrf,  and  so  is  co(V). 

Claim: 

co(V)  n  (R+)rf  =  0  (59) 

where  Rq  is  the  set  of  nonnegative  real  numbers. 

In  fact,  if  (59)  does  not  hold,  then  we  can  find  a  point  y  in  the  boundary  of  co(V),  which  falls 
inside  (Mj )d .  Invoking  Lemma  10,  y  is  a  convex  combination  of  d  points  in  V.  Therefore  there 
exist  (Ji,.  01 , . . . ,  0d,  vl , . . . ,  vd  such  that 

d  d 

y  =  '52akK,„,e*(vk)  6  (Ro)d,  ak>0,  Y,ak=zl-  (60) 

k— 1  k= 1 

Recalling  (56),  (57),  then  (60)  implies  that  (48)  holds  for  the  corresponding  G(wk,0k)  and  vk  = 
y/ OkVk .  This  violates  (6)  by  Lemma  11,  so  the  claim  is  proved. 
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co(V)  and  (Rq  )d  are  disjoint  closed  convex  sets  in  Rd,  and  co(V)  is  compact.  Therefore,  by 
Lemma  9  there  exists  x  =  (xi, . . xj)  G  Rd,  a,/3  G  R  such  that 

(x,A)  <  a  <  f3  <  (x,y)  VA  G  V,  y  G  (R^ )d  (61) 

Given  the  special  structure  of  the  cone  (Rq  )d,  we  can  choose  (3  =  0  and  x\  >  0, . . . ,  xj  >  0.  Using 
(57),  we  have 

xial,e(v)  +  •  •  •  +  xd°ifi(v)  <  «  <  0  Vo;  G  [-x, x],  6  G  B©,  v  G  Cm,  ||n||  =  1  (62) 

Let  X  =  diag[x\Imx  . .  .Xdlmd\  which  belongs  to  then  (56)  and  (62)  give 

\\X2G(u,6)v\\2  -  \\X*v\\2  <  a  <  0  Vw,  6,  v  (63) 

It  follows  that  7  :=  max^d  (^X^G(u>,  6)X~^  <  1.  For  any  LTI  9(ejM)  taking  values  in  B0,  and 
any  LTV  $  in  B$,  we  obtain 

XH(M  *  0)X~2  <  X2G(u,e(ej“))X~2  <7  (64) 

OO 

This  implies 

(i  -  XH(M  *e)X~12Y1  <— !—  (65) 

v  /  1  —  7 

which  leads  to  a  uniform  bound  on  ||  (/  -  $(M  *  0))_1  j|.  This,  together  with  (55)  and  assumption 
(20)  gives  a  uniform  bound  on  ||(J  -  AM)-1||. 

□ 

7.1  Proof  of  Proposition  5 

It  remains  to  prove  that  (29)  implies  (30).  (29)  implies  that 

max  infd(XM(i^,...,wASr-1)  <  1  (66) 
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For  fixed  <uq  , . . . ,  u>d,  clearly  there  exists  e(uq , . . . ,  uj)  >  0  such  that 

inf  afXMfa , . . . , Wd)!-1)  <  1  (67) 

\l>X>cI 

Since  the  oq,  . .  .,a >d  vary  in  a  compact  set,  and  M  is  continuous,  it  follows  that  a  fixed  e  can  be 
found  satisfying  (67)  across  oq,. .  .,o><£.  Now  assume  that  (30)  does  not  hold.  This  means  that  the 
LMI  condition  across  frequency 

X(u)  >  0,  M*(lo)X(cv)M(u)  -  X(u)  <  0  (68) 

does  not  have  a  solution  X  £  Xm,  i.e.  it  cannot  be  satisfied  with  the  X®  part  constant  across 
frequency.  As  a  consequence,  the  family  of  sets  over  u 

X®  :=  e  X*  :  3X  =  diag[X*,  X8(w)},  el  <  X  <  M(ejuj)*X M(ej“)  -  X  <  oj  (69) 

have  empty  intersection  for  u  ranging  in  [— 7r,7r].  These  sets  are  convex  and  compact,  and  by 
normalizing  the  last  block  of  X®  to  /,  they  are  in  a  d  —  1  dimensional  space.  From  Lemma  8,  there 
exist  d  such  sets  with  empty  intersection. 

Therefore  there  exist  d  frequencies  uq .,Ud  such  that  a  [X M(ju>k)X~l)  <  1  cannot  be  sat¬ 
isfied  with  el  <  X  <  f  and  a  common  X®.  This  contradicts  the  fact  that  e  satisfies  (67)  for  all 
frequencies  uq, . . .  ,0Jd- 


□ 


7.2  Proof  of  Theorem  6 

The  sufficiency  of  condition  (30)  for  uniform  robust  stability  over  J~{v)  can  be  proved  with  very 
minor  modifications  to  the  proof  given  in  [14];  the  details  are  omitted.  We  now  outline  the  necessity 
proof,  which  is  essentially  a  combination  of  the  methods  applied  in  [14]  to  deal  separately  with  slowly 
varying  perturbations  and  “arbitrarily  fast”  perturbations. 
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Necessity:.  For  simplicity  we  assume  the  structures  0  and  $  consist  only  of  full  blocks, 
€  £(1™‘),  6  We  have  A  =  diag[$\, . . .,  $d,9d+\,  ■  ■  • ,  ^+f]- 

Assume  now  that  the  system  does  not  satisfy  (30);  we  will  show  uniform  robust  stability  is 
violated  for  every  v  >  0. 

Step  1:  We  first  use  Proposition  5  to  show  the  existence  of  d  frequencies  uq, . . .  ,u>d  such  that 

(70) 

At  this  point,  we  will  perform  a  second  finite  augmentation  of  the  problem  in  the  style  of  Section  3. 
For  notational  simplicity  a  matrix  form  will  not  be  written,  but  it  follows  from  [2]  (or  also  from  the 
proof  in  Section  3)  that  there  exists  a  finite  horizon  N  and  vectors  v1, . .  .,vN  6  Cdm+p)xd  which 
“destabilize”  an  augmentation  of  the  matrix  M  in  the  following  sense. 

For  each  n  £  1 . . .  JV,  partition  the  vector  vn  as  ( vk)n  in  accordance  to  the  blocks  of  M;  let 
M{e^k){vk)n  =  ( yk)n ;  the  following  inequalities  hold: 

£!|((/n|faE||(<'‘)")I||  '  =  !••.<(  pi) 

k,n  k,n 

E||((»‘)”wf  2£||((»‘r)«/f  { Izlz*  (72) 

n  n  ^ 

Note  that  all  “energy  can  be  shared”  between  the  repetitions  n  of  this  second  augmentation.  As  for 
the  first  augmentation  in  k,  only  the  $  blocks  “transfer  energy”  as  before.  Heuristically,  this  second 
augmentation  will  be  used  with  frequencies  close  to  the  Wk  s.  The  transfer  of  energy  between  these 
frequencies  can  be  achieved  by  slowly  varying  operators. 

Step  2:  Fix  u  >  0  and  e  >  0.  For  6  to  be  chosen  below,  construct 

N-l  N- 1 

vk(t)  =  ^(ufc)"ei^+2^15),  yk(t)=  J2(yk)nej{u'k+a*L6)  (73) 

n=0  n= 0 

By  continuity  of  M(e^)  around  Uk,  yk  =  Mvk  +  ek  +  wk,  where  wk  is  a  transient  term  in  I2, 
and  ek  is  a  sum  of  sinusoids  such  that  its  power  ||  •  || -p  satisfies  \\ek\\ p  <  e||ufc||^>,  for  small  enough  6. 
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Also  choose  6  small  enough  so  that  2  sin(  |)  <  v ,  and  so  that  the  intervals  [wk,  Wk  +  £]  are  pairwise 
disjoint. 

From  (72)  and  Lemma  2,  for  each  /,  k  we  can  find  0k+f  £  F{y)  such  that  Ok+jyd+,  -  vd+r 
Since  the  Fourier  supports  of  ( yk,vk )  are  disjoint  for  different  k ,  with  the  same  methods  as  in 
Lemma  2  it  can  be  shown  that  there  is  a  single  9d+f  G  F{v),  satisfying 

Od+fVd+f  =  vd+f,  k  =  l...d  (74) 

Now  define  v(t)  =  J2k=i  vO-)  =  ^2k= l  Vk(t)-  We  have  y  =  Mv  +  e  +  w,  w  £  h,  ||e|| j,  = 

0(e)  ||u||-p,  and  the  definition  of  the  vk ,  yk,  implies  from  (71)  that 

INI2  >  INI2,  l=l...d  (75) 

This  implies  we  can  construct  contractive  LTV  operators  l  =  1 . .  .d  (with  no  prescribed  rate 
of  variation,  since  v,  y  are  not  confined  to  a  “thin”  frequency  band)  such  that  §iyi  =  v;.  Setting 
A  =  diag[$i,  &d+ 1 ,  •  •  • ,  Od+F\  gives,  using  (74),  Ay  =  v.  Therefore 

||(i"  —  AM)v\\-p  =  HAeUp  =  0(e)  \\v\\v  (76) 

Since  e  is  arbitrary,  uniform  robust  /2-stability  is  violated  within  the  class  A  =  diag[$,Q],  0  G 
T(v).  Since  u  is  arbitrary,  this  completes  the  proof. 

□ 
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